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O Abstract 
O ' 
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i-C . We construct a dilatonic two-dimensional model of a charged black hole. 

<D 

£Li - 

The classical solution is a static charged black hole, characterized by two 
parameters, m and q, representing the black hole's mass and charge. Then we 
J> ■ study the semiclassical effects, and calculate the evaporation rate of both m 

g ■ and q, function of these two quantities. Analyzing this dynamical system, 

■ we find two qualitatively different regimes, depending on the electromagnetic 

coupling constant qa- If the latter is greater than a certain critical value, the 
O" 1 charge-to-mass ratio decays to zero upon evaporation. On the other hand, 

bJ[). for gA smaller than the critical value, the charge-to-mass ratio approaches a 



X 



non-zero constant that depends on gA but not on the initial values of m and 



I. INTRODUCTION 

Dilatonic two-dimensional models are very useful in studying various aspects of black 
holes (BHs). Callan, Giddings, Harvey and Strominger (CGHS) |]J introduced a dilatonic 
two-dimensional model to investigate the evaporation of a BH. They further used the two- 
dimensional model to explore the final stage of the evaporation process. The latter, more 
ambitious, goal of clarifying the endpoint of evaporation has proved to be illusive and hard 
to achieve ||. Nevertheless, the two-dimensional model provided a useful description of the 
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evaporation process in its initial phase, i.e., as long as the BH's mass is sufficiently large. 
Since then, many authors studied various aspects and various variants of CGHS's model. 

The two-dimensional model studied by CGHS was that of an uncharged BH. The main 
objective of this paper is to develop an analogous model to describe the evaporation of an 
electrically charged two-dimensional BH. Our main motivation for studying charged BHs 
stems from the basic features of the classical, macroscopic, 4-dimensional BH solutions. The 
spacelike singularity of the Schwarzschild geometry disappears when an electric charge is 
added to the BH, and an inner horizon (IH) forms instead. Remarkably, the same situation 
occurs when an angular momentum is added to the BH (instead of an electric charge): 
In the Kerr solution, too, there is an inner horizon and no spacelike singularity. This 
observation is very relevant to reality, because realistic astrophysical BHs are believed to be 
rapidly spinning || . Spherically symmetric charged BHs thus provide a useful toy model for 
exploring various aspects related to the inner structure of the more realistic, spinning BHs 
(see e.g. Ref. §). 

In the last decade several authors investigated two-dimensional models of charged BHs 
(see e.g. [|], ||, J7|, || ). McGuigan, Nappi, and Yost |5j studied such a classical model, 
with a dilaton coupling to the electric field. They considered a coupling e -2 * of the dilaton 
field to the electromagnetic term in the action (see below). They found that, just as 
in the four-dimensional case, charged BHs admit an inner horizon instead of a spacelike 
singularity. However, their model does not include semiclassical effects, which are necessary 
for describing the BH's evaporation. 

Later, Nojiri and Oda (NO) considered a slightly modified dilatonic model, in which 
the dilaton coupling to the electric field is e 2< ^. In their model there is a large number N 
of chiral fermion fields (instead of the iV scalar matter fields considered in Ref. ) , which 
couple to both the curvature and the electromagnetic field. NO first studied the structure 
of the classical BH solution. Then, generalizing the method used by CGHS to the charged 
case, they added two effective semiclassical correction terms to the action, and derived from 
them the field equations at the semiclassical level. They have been able to solve some of the 
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semiclassical field equations, which allowed them to analyze various aspects of the resultant 
semiclassical charged BH solutions. 

Nojiri and Oda considered the dilaton coupling e 2 ^ rather than e _2< ^, primarily because 
the former coupling makes some of the equations easier to solve. As it turns out, however, 
there is a remarkable difference between the two models, already at the classical level. The 
charged BHs with the coupling e 2< ^ do not admit an inner horizon || (instead they usually 
admit a spacelike singularity, just as in the uncharged case). Since our main motivation for 
considering charged BHs is to mimic the inner structure of the four-dimensional spinning 
BHs, we find it important to elaborate on the charged model with the dilaton coupling e~ 2 ^. 
(Note also that e~ 2 ^ is the coupling which emerges as the effective action in the low-energy 
limit of string theory. ||) 

Motivated by these considerations, in this paper we start from the model developed 
by NO, and modify the dilaton coupling to the electromagnetic field from e 2 * to e~ 2<? \ 
The resultant field equations are harder to solve at the semiclassical level. Nevertheless, 
it is possible to solve the equations describing the semiclassical effects in the adiabatic 
approximation, i.e. in the approximation where the background geometry is described by 
the static, non-evaporating, BH. [[5] This approximation appears to be valid as long as the 
black hole is macroscopic, in which case the evaporation rate is small (namely, the relative 
change in the mass or charge during a dynamical time scale is << 1). It is this macroscopic 
domain which will concern us throughout this paper. Thus, on the background of a classical 
static BH solution (with given mass and charge), we solve the semiclassical equations and 
derive from them the semiclassical fluxes of both energy-momentum and charge. This allows 
us to determine the evaporation rate of both the mass and charge of the BH. 

We start in section II by writing the classical action and the corresponding field equations. 
We define new variables, which are used to simplify the equations and to present their general 
solution. This general solution is a two-parameter family of static black hole solutions, 
uniquely characterized by the two parameters m and q, which represent the BH's mass and 
charge. This classical solution was already given in Ref. H, though in different coordinates. 
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Here we construct the classical solution in double-null coordinates, which are more suitable 
for the subsequent semiclassical calculations. We also extend the classical solution to include 
an outgoing (or ingoing) null fluid. (This extension is useful for describing the geometry of 
the evaporating BH at large distance from the horizon.) 

In section III we consider the semiclassical effects. Following NO, we analyze the semi- 
classical effects by adding two effective correction terms to the classical action, expressed 
in terms of two new variables Z and Y (these variables describe the semiclassical fluxes of 
energy-momentum and charge, respectively). From this action we derive the semiclassical 
field equations. Then we solve the field equations for the two quantum variables Z and Y, 
assuming a background geometry of a static, classical BH. This solution yields an explicit 
expression for the fluxes of energy-momentum and charge, at any location (both outside 
and inside the BH), as a function of m and q. Based on these fluxes, in section IV we 
calculate the rate of evaporation of the mass and charge. We obtain a closed system of two 
first-order equations, describing m and q as functions of m and q, where an overdot denotes 
a derivative with respect to the external time. Both m and q are found to be negative, 
as one should expect. We then analyze this dynamical system, and find two qualitatively 
different regimes. If the electromagnetic coupling constant g A is larger than a certain critical 
value, the charge-to-mass ratio decays to zero upon evaporation (as a certain power of the 
mass, which depends on g A ). On the other hand, for g A smaller than the critical value, the 
charge-to- mass ratio approaches a non-zero constant, that depends on g A but not on the 
initial values of m and q. 

In section V we summarize and discuss our results. It should be emphasized that no 
attempt is made in this paper to investigated the final stage of evaporation. The goal here is 
to study the semiclassical evolution of the BH in the macroscopic phase, i.e. as long as the 
BH's mass is much larger than a certain mass. It is this phase in which the above mentioned 
adiabatic approximation is valid. In section V we further discuss this validity criterion for 
the adiabatic approximation, and find the range of mass values for which this approximation 
can be used. 
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II. CLASSICAL SOLUTIONS 



A. field equations 

We start with the classical action S c given by NO |J, which includes an electromagnetic 
field coupled to charged matter represented by N left-handed chiral fermions: 



S c = — [ d 2 x ^{e~ 2 *(R + 4(V0) 2 + 4A 2 
2% J 



N 



--jP-Yt&ffiPr-iA^,}. (1) 

3=1 

Here 6 is a dilaton field, = (*>) are the N left-handed chiral fermion fields, qa is the 
electromagnetic coupling constant, and denotes a covariant derivative. The Maxwell 
tensor F^ u is given by F^ v = d^A u — d u A^, with F 2 = F^ V F^ V . The coupling of the dilaton 
to the electromagnetic field term is written here, in a quite general form, as e a< ^. In Ref. ||, 
NO only analyzed the case a = 2. Here, for the reasons explained in the Introduction, we 
shall consider the case a = —2. 

Following NO, we use the light-cone gauge for the electromagnetic field, namely 

A u = . (2) 

We also use double-null coordinates u, v, with 

Quv 2^ ^ i 9uu Qw 0. (3) 
The action (P (with a = — 2) then reduces to 

Sc=i LJ d2x {e ~ 2 ^ 4/) <™ ~ & ^> v + 2avp) 

A ' N 

+V 2( ^^ + |E^i}' ( 4 ) 

1 3=1 

where F uv = d u A v . The Einstein equations are given by 
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= T VV = e- 2 \A P ^^ v - 2<f> iV0 ) + - Y^WjWi ~ WW) 

3=1 

1 N 

+-A v j2rw, (5) 

3=1 

= T uu = e-^(4 PiU iU - 2^ uu ) , (6) 



= T uv = e- 2<f> (2^ uv - 40, u 0,„ - AV) + A e -2(0 + p) F ^ (7) 



The dilaton equation of motion is 



= -40,™ + 40>,„ + 2p iUV + X 2 e 2p + ^e~ 2p F 2 v , (8) 

9a 



the Maxwell equations are 

8 i N 

= \d v {e~ 2 ^F uv ) + -YVV; , 
^ i^i 

= ^^(e- 2 ^F u „) , (9) 
9a 

and the fermion fields satisfy 

= d u iPj . (10) 

We shall consider here solutions free of any classical matter. [] That is, we shall only 
consider here the trivial solution 

lfc = (11) 
to Eq. (0). The Maxwell equations then reduce to 

= d u (e~ 2 ^F uv ) = d v {e~ 2 ^F uv ) , (12) 

namely, 



1 Presumably the charged BH was created by the collapse of the fermionic matter. Here, however, 
we are interested in the evaporation of the BH and not in its creation. 
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(V2/g A ) e~ 2 ^ +p) F uv = const = Xq. (13) 

Substituting Eqs. (|T0| ) and (|l3j) in the above system([5]-|D, we obtain a simpler closed system, 
which includes the three Einstein equations 

= T uv = e- 2 *(20, ra - 4<^0,„ - A 2 e 2 ' + A 2 g 2 e 4 * +2 ") , (14) 

= T vv = e- 2 *(4p>^ - 2<f> >vv ) , (15) 

— T uu = e- 2 *(4p in jU - 20, nu ) , (16) 

and the dilaton equation 

= -40,™ + 40>,„ + 2p^ + A 2 e 2 " + A 2 g 2 e 4 * +2 " . (17) 

Note that the two equations ([14]) and (|I~7D - to which we shall refer as the evolution equations 
- are hyperbolic, and are hence sufficient for determining the evolution of the two unknowns 



p and (f> from prescribed initial data. [The other two equations (|T5| ,[Tj]) - the constraint 
equation - are consistent with the evolution equations: Any solution of the evolution equa- 
tions whose initial data are consistent with the constraint equations, will also satisfy the 
constraints in the entire domain of dependence.] 

B. New variables 

To further simplify the analysis, we define the new variables 

R = e - 2 * , S = 2(>-0). (18) 
The evolution equations then reduce to 

Ruv = X 2 (q 2 /R 2 - 1) e 5 , (19) 

S, uv = (-2X 2 q 2 /R 3 ) e s . (20) 
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The two constraint equations also take a simple form: 

= T uu = R :UU R U S^ U , (21) 
= T vv = R^ vv R V S^ V . (22) 



The two evolution equations (|T^j20| ) can be expressed in a compact form as 

R :UV = F(R)e s , S :UV = F R e s . (23) 

where 

F(R) = X 2 (q 2 /R 2 - 1) . (24) 

As it turns out, the system (B^) is rather universal, as various general-relativistic models 
(e.g. several two-dimensional BH models, the three-dimensional BTZ model, and various 
spherically-symmetric four- dimensional models) satisfy the same form of hyperbolic system, 
with each model having its own function F(R). For example, in the spherically-symmetric 
four- dimensional model of a charged BH with (or without) a cosmological constant, if one de- 
fines R = r 2 and e s = rg uv (where r is the area coordinate and u, v are two null coordinates), 
the two Einstein evolution equations take exactly the form ( |23|) with 



F(R) = aR 1/2 + bR- 1/2 + cR~ 3/2 (four - dimensional) , (25) 

where a, b, c are constants (a and c represent the contributions of the cosmological constant 
and charge, respectively). 

The non-linear hyperbolic system (^) [for a rather generic function F(R)} may serve as 
a useful mathematical toy model for studying various aspects of the theory of black holes, 
like gravitational collapse, singularity formation, and the no-hair principle. This, however, 
is beyond the scope of the present paper. Here we shall merely use the R, S variables to 



simplify the analysis, as the system ( p3|) does not include first-order derivatives. We shall 
also use a few general features of this system - e.g. the form of its static black-hole solutions. 
The generic solution of Eq. (|2"3"D does not necessarily satisfy the constraint equations 



(^T|,^). Of course, we are primarily interested here in the subclass of solutions which do 



satisfy the constraint equations, to which we shall refer as the vacuum-like solutions. Apart 
from the gauge freedom (i.e. the freedom to re-parametrize each of the two null coordinates 
u,v), this subclass is a one-parameter family of solutions [for a given F(R) ], parametrized 
by the mass. These vacuum-like solutions are, in fact, the static black- hole solutions of 
the model. (For example, in the four- dimensional spherical electrovacuum case, these are 
the RN-deSittre family of solutions.) In the context of the specific model considered in 
this paper [with F(R) of Eq. fl24]) 1, the vacuum-like solutions are nothing but the two- 
dimensional static electrovacuum solutions. We shall construct these static solutions, in 
double-null coordinates, in the next subsection. 

One may also be interested in the wider class of solutions to the evolution equations 



(p3|), which do not necessarily satisfy the constraint equations (|2T| , |22D . Such solutions may 
be interpreted as spacetimes perturbed by ingoing and/or outgoing null fluids, leading to 
non- vanishing contributions to T vv and/or T uu . We shall name such solutions as radiative 
solutions. For example, in the context of spherically-symmetric, four-dimensional, charged 
BHs, the mass-inflation solutions introduced in Ref. || belong to this class of radiative 
solutions [with the function F(r) of Eq. (P5|)]. Note that although a radiative solution does 
not satisfy all the vacuum field equations, its evolution from Cauchy or characteristic initial 
data is completely determined from the (vacuum!) evolution equations, which form a closed 
hyperbolic system. 

An important subclass of the radiative solutions are those which satisfy one of the con- 
straint equations, but not the other one. Such solutions may be interpreted as spacetimes 
with a null fluid flowing in either the outgoing or ingoing direction. We shall refer to these 
solutions as the Vaidya-like solutions. For example, the geometry in the weak-filed region 
(i.e. at large R) of an evaporating BH can be well approximated by an outgoing Vaidya-like 
solution. In Appendix A we describe the construction of the general Vaidya-like solution in 
double-null coordinates. 
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C. The static black-hole solution 



The general solution of the above system |22J) (i.e. both the evolution and constraint 
equations) is a family of two-dimensional static, RN-like, black-hole solutions uniquely char- 
acterized by their mass and charge. This general solution was presented in Ref. using 
Schwarzschild-like coordinates. For the analysis below we shall need the solution in double- 
null coordinates. We shall first describe the construction of this solution for a general 
function F(R), and then restrict attention to our specific model, i.e. F(R) = X 2 (q 2 /R 2 — 1). 

For a general function F(R), we define 

H{R) = - F(R')dR' . (26) 

The static, vacuum-like, solution only depends on the spatial coordinate, which we denote 
x. We choose an Eddington-like gauge, such that x = v — u. The solution is then given by 

e s = H(R) , R iX = H(R) , x = v — u (H > 0). (27) 

From the above definitions of R and S, the metric function g uv is given by 

-2g uv = e 2p = H/R. (28) 

Note that this Eddington-like solution is only valid in the region outside the BH where 
H(R) > - this is the region which will primarily concern us in this paper. The solution 
exhibits a coordinate singularity whenever H vanishes [where g uv vanishes, and so does det(g) 
]. The lines H = correspond to the horizons of the BH. These include the event horizon 
(EH), and [for functions H(R) which admit more than one zero] also the inner horizon 
and/or the cosmological horizon. In the region inside the BH where H(R) is negative, the 
solution can be expressed in a very similar form - see section III. Note that the solution 
includes one free parameter - the integration constant in Eq. (^) - which is related to the 
BH's mass (see below). 

In the specific model considered in this paper, for which F(R) = X 2 (q 2 /R 2 — 1), we write 
H(R) in the form 
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H{R) = X 2 (R - 2m + q 2 /R) (29) 

(for notational convenience we take here the integration constant to be —2mA 2 ). The 
vacuum-like solution is thus uniquely determined by the two parameters m and q, which are 
related to the black-hole's mass and charge, respectively. The root structure of the function 
H(R) depends on the ratio between q and m. In this paper we shall consider non-extreme 
charged black-hole solutions, i.e. solutions with m > q > (the restriction to positive rather 
than negative q does not cause any loss of generality). The equation H(R) = then has 



two roots, at R± = m ± \Jm? — q 2 , where R + corresponds to the EH and corresponds to 
the IH. The function H(R) is positive outside the BH, i.e. at R > R + (and also at R < 
but this range will not concern us in this paper), and negative between the two horizons. 
As was mentioned above, the solution ( p7|) is only valid outside the BH, and a similar one, 
valid inside the BH, is given in section III. 

The surface gravity k + of the EH is defined by 

k + = \{H, R ) R+ = A 2 (l - q 2 /R 2 + )/2 . (30) 

For later convenience, we also express k + in other useful forms: 

k+ = A 2 (l - m/R + ) = A 2 [(m 2 - q 2 ) l/2 /R+] . (31) 

It is remarkable that in the two-dimensional case, unlike the situation in four- dimensional 
BHs, n + only depends on q/m, and not on the BH's size. One explicitly finds 

-r 



1- 1 + y/l-(q/m)' 



(32) 



Note that k + is a decreasing function of q/m, and for < q/m < 1 it takes the values 
< k + < A 2 /2. 

From Eq. (|28|) , the metric function g uv is given by 

-2g uv = \ 2 {l-2m/R + q 2 /R 2 ). (33) 

At large R, this becomes —2g uv = A 2 . It is useful to introduce new null coordinates 
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u = Aw, v = Xv , 



(34) 



such that 



1 - 2m/ R + q 2 /R 2 , 



(35) 



which yields the desired asymptotic behavior, —2gnv = 1, at large R (this is also the type of 
gauge used by CGHS in the uncharged case). We shall refer to the u,v coordinates as the 
asymptotically-flat double-null coordinates. 

Since R = e~ 2r ^ is dimensionless, the two parameters m and q are dimensionless too. 
These two parameters are proportional to the physical mass and charge of the BH, and 
we may refer to them as the dimensionless mass and charge, respectively. [One of the 
motivations for this association is the similarity of the metric function g^a in Eq. (p5|) to 
its counterpart in the standard four-dimensional Reissner-Nordstrom solution.] Throughout 
this paper we shall often refer to m and q simply as the BH's mass and charge (with some 
abuse of the standard terminology). To avoid confusion, we shall denote the dimensionful, 
physical, mass and charge of the BH by M and Q, respectively. The physical mass M, 
which is the total energy content of the system, is encoded in the asymptotic behavior of the 
metric tensor at large distance. Since this large distance corresponds to the limit R — > oo, 
the physical mass will not be affected by the term q 2 /R 2 in the metric (^). Hence M 
must be a function of m and A solely. We can deduce M(m, A) from the form of the mass 
parameter M in the uncharged case, studied by CGHS. Comparing the static solutions of 
the two models, one finds 



2 To relate M in the CHGS model to our notation, we can compare the value of the dilaton field 
at the EH. In our notation we find at the EH (for q = ) e~ 2 ^ = R = 2m, whereas in CGHS we 
find at the EH e" 2 ^ = M/A; cf. Eq. (11) therein. 



M = 2Xm . 



(36) 
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The relation between q and Q may be revealed by comparing our static solution to that 
given (in different coordinates) in Ref. ||. This comparison shows that Q is proportional to 
Xq , in analogy with Eq. (0). We have not clarified the constant relating Q to Xq (however, 
the explicit expression for Q is not required for the analysis below). 

III. SEMICLASSICAL CORRECTIONS 

A. Semiclassical field equations 

In order to study the evaporation of the BH, we must consider the semiclassical contribu- 
tion to the energy- momentum tensor and to the electromagnetic current. CGHS |l| showed 
that the semiclassical contributions can be treated by adding an effective correction term to 
the classical action. Nojiri and Oda |J extended this method to the charged case (they also 
modified the action by considering fermionic rather than scalar matter fields). They found 
that the semiclassical contributions coming from the conformal anomaly can be represented 
by adding two correction terms to the classical action S c : 

S = S c + S p + S x , (37) 

where 

S P = ~ J #x^{--{VZ)* + yj-ZR} (38) 

and 

S x = i J d 2 x {-iv^(VF) 2 + ^Ey^F^} . (39) 

The term S p , which emerges from the trace anomaly, contributes to the energy-momentum 
tensor. The second term S x comes from the chiral anomaly, and it contributes to both the 
electromagnetic current and the energy-momentum. Note the presence of two new variables, 
Z and Y, in the correction terms. These two degrees of freedom were introduced by NO in 
order to allow the representation of the semiclassical effects by local correction terms 0. 
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Expressing the metric in the double-null form (|3]) and using the electromagnetic gauge 
(0), the modified Einstein equations become 

= T VV = e- 2lf, (4p, v( J), v - 2<P tVV ) + f vv , (40) 

= T uu = e- 2 ^(4p in iU - 20, nu ) + f uu , (41) 

= T uv = e-^(20, ra - 40, u )t) - X 2 e 2p + q 2 e wp ) + f uv . (42) 

Here, T uu , T vv , and T uv represent the semiclassical contributions to T uu , T vv , and T uv , 
respectively, which are given by 

%v = \z% - \[^P,vZ >v + l\[^Z, vv + l -Y 2 , (43) 

Tuu = \z% - ^P,uZ, u + \\j^Z, uu + l -Y 2 , (44) 

T U v V 48^'™ ' (^^) 
The Maxwell equations also get semiclassical source terms: 

= \d u {e~ 2 ^F uv ) - V2Nd u Y , (46) 
9a 

= \d v (e~ 2 ^F uv ) - V2Nd v Y . (47) 
9a 

The dilaton and matter equations of motion (||) and (|10"D are not modified. [As before, we 
consider the vacuum solution ( |Tl~| ) to Eq. ([10]).] The variables Z and Y satisfy the field 
equations 

= -2Z iUV + \j^P,uv , (48) 

= -2Y UV - V2NF UV . (49) 
Motivated by the classical relation flT3|), we define q by 



\q = (V2/g A )e- 2 ^ +p) F uv . (50) 
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Note that q is no longer constant: It evolves due to the semiclassical source terms in the 



modified Maxwell equations (|46| , |4"T| ). These two equations can immediately be integrated: 



q = q + (g A /4\)VNY 



(51) 



where go is an arbitrary integration constant. The equation for Z can be integrated too: 



N 




Z = \l —p + Z v {v) + Z u (u) 



(52) 



where Z v (v) and Z u (u) are arbitrary initial functions. 

To simplify the equations, we define K = iV/24, g = V3gA , and rescale Y and Z as 



z = Z/V2K , y = Y/V2K. 



(53) 



Equations ( f49| , |5"T| , |52"D then become, respectively, 



jj.l 



-Vqf u 



(54) 



q = q + {Kg/X)y. 



(55) 



and 



z = p + z v (v) + z u (u) 



(56) 



where again z v (v) and z u (u) are arbitrary initial functions. The semiclassical contributions 
to the stress-energy now read 



T uu K {^z u 2p u z <u -\- z uu -|- U tU j 



T 



-Kp, 



Using Eq. (|50D, we can rewrite the field equation for y in terms of q: 



y, uv = -Xg e 2( ^g . 



(57) 
(58) 
(59) 



(60) 
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The closed system of semiclassical field equations is composed of Eq. ( [LTD for the dilaton, 
the three Einstein equations ( [TO[]42| ), as well as Eqs. ( |55| - |6"0D . 

Finally, we write the system of semiclassical field equations in the R, S variables. The 
dilaton equation ( |TTD and the Einstein equation fl4~2| ) [with Eq. (|59|)1 yield 



Ruv = X 2 (q 2 /R 2 - 1) e s -Kp tV 



(61) 



S, 



(2X 2 q 2 /R 3 ) e s + Kp, uv /R 



(62) 



where p = (S — logi?)/2. The semiclassically-corrected constraint equations take the form 



R,uu R,uS^u T uu o 
R,w R.vS d ~\~ T vv 



(63) 
(64) 



Again, this system is supplemented by Eqs. (|55] - [K]). Equation ( |60| ) can be re-expressed 
using the R, S variables as 



y,uv = -Xg (e s /R 2 ) q ■ 



(65) 



B. Semiclassical fluxes outside the black hole 

We turn now to analyze the evolution of the quantum variables Y and Z, in order to 
obtain the semiclassical fluxes. To that end we use the adiabatic approximation. Namely, 
we view Y and Z as test fields living on the background described by the static classical 
solution (with fixed m and q). 

We first calculate the semiclassical electric currents outside the BH. Using Eq. (|27|), we 
rewrite Eq. ( |63D as 

V,uv = -XgqH/R 2 . (66) 
We now integrate this equation with respect to v, recalling dv = dx = dR/H: 
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V,u = -Xgq J (H/R 2 )dv = -Xgq J R~ 2 dR = Xgq/R + J u {u) . (67) 

Similarly, we find for y v (recalling du = —dx): 

V v = -Xgq J (H/R 2 )du = Xgq J R~ 2 dR = -Xgq/R + J v {v) . (68) 

The integration constants, i.e. the functions J u {u) and J v (v), are to be determined from 
the initial conditions. Since we assume no ingoing current is coming from past null infinity, 
we must set J v (v) = 0. Also, regularity at the EH requires that y >u vanishes there, which 
implies J u {u) = —Xgq/R + . Therefore, 

y !U = Xgq{l/R-l/R + ) , y, v = -Xgq/R, (69) 

and from Eq. (|55|) we obtain 

q iU = Kg 2 q(l/R-l/R + ) , q, v = -Kg 2 q/R. (70) 

The fluxes T vv and T uu , Eqs. (|57|j58|), can be expressed explicitly by means of Eq. (p6|): 



f uu = K[(p tUU - p 2 u ) + y 2 u + z u (u)} , (71) 
f vv = K[{p tVV - p 2 ) + y 2 + z v (v)} , (72) 

where z v (v) = z v>vv + z 2 v and z u {u) = z UyUU + z 2 u . Substituting Eq. (p5|), we find 

f uu = K[(p, uu - p%) + X 2 g 2 q 2 (l/R - l/R + ) 2 + z u (u)] , (73) 
f w = K[(p, vv - p 2 v ) + X 2 g 2 q 2 /R 2 + z v (v)\ . (74) 

The functions z v (v) and z u (u) are to be chosen such that no influx is coming from past null 
infinity, and the outflow is regular at R = R + , that is, 

f vv (R = oo) = , f uu (R = R+) = 0. (75) 

In the static classical background we have p = (1/2) log(H/R), so 

- P,u = p, v = P, x = Hp, R = (R/2)(H/R) iR (76) 
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and 

P,uu = P,w = P, xx = H[{R/2){H/R), R ) >R (77) 

Note that p tUU and p jVV (as well as p tUV ) vanish both at R + and at R — oo. On the other 
hand, p >u and p tV vanish at R — oo, but at R = R + they get a finite value, 

- P,u = P,v = «+ , (-R = (78) 
In order for T„„ to vanish at R = oo, we choose z v = and obtain 

T TO = #[(P, W - P 2 V ) + (Xgq/R) 2 ] • (79) 
Also, the demand that T uu vanishes at R = R + yields z u = k 2 + , namely, 

f uu = K[k\ + (p tUU - pi) + [Xgq(l/R - 1/R+)} 2 ] . (80) 

C. Semiclassical fluxes inside the black hole 

Before we discuss semiclassical effects inside the BH, we need to extend our classical 
solution for the static black-hole background to the internal region. Clearly, Eq. Q27D as it 
stands is not valid at R_ < R < R + , where H is negative. The internal solution in double 
null, Eddington-like coordinates is obtained form Eq. (|27]) by minor changes of sign. The 
main difference is that, inside the BH the variable x (the only variable on which the solution 
depends) is temporal rather than spatial, namely, x — v + u. The internal solution takes 
the form 

e s = —H(R) , R >x = H(R) , x = v + u (R_ < R < R+) . (81) 

Correspondingly, the metric function g uv is given by 

-2g uv = e 2p = -H/R. 

We can now repeat the calculations of the semiclassical fluxes. The initial conditions for 
the outgoing fluxes at v — > — oo are the same as in the external problem: Both the energy 
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and charge outfluxes must vanish at R = R+, in order to achieve regularity at the EH. The 
initial conditions for the ingoing fluxes at u —>■ — oo (the EH) are dictated by continuity: 
At the EH, both q )V and T vv must continuously match the corresponding quantities in the 
external region R > R + (recall that the coordinate v is regular at the EH). 

The calculation now proceeds in a way completely analogous to the external semiclassical 
calculations of the previous subsection, except for a few changes of sign. For example, 
when calculating the charge fluxes, one must recall that e 5 = —H and du = dx, and as a 
consequence q }U changes sign (but not q tV ): 

q :U = -Kg 2 q(l/R-l/R + ) , q, v = -Kg 2 q/R (R_ < R < R + ) . (82) 

The energy-momentum fluxes are given by Eqs. (|79|) and fl8~0|) without any change (recall, 
though, that now p is given by e 2p = —H/R). 

IV. EVAPORATION OF THE BLACK HOLE 

A. Evolution of m and q 

In order to calculate the rate of change of m and q, as measured by an observer at future 
null infinity (FNI), we need to evaluate the outgoing fluxes at the limit R — > oo. For brevity 
we denote the u -derivatives of m and q at FNI by an overdot. At this limit Eq. (|70|) reads 

q = -Kg 2 q/R + . (83) 

Taking the large- R limit in Eq. (|80|), we find 

f uu = K[k 2 + + (Xgq/R + ) 2 } = KX 2 [X 2 (m 2 - q 2 ) + g 2 q 2 ]/R 2 + (FNI) . (84) 

The relation between rh and T uu is most easily expressed in terms of the BH's physical mass 
M = 2Xm and the asymptotically-flat null coordinate u = Xu (see section II): The change 
in the Bondi mass M is simply (minus) the integral of T^u with respect to u along FNI; that 
is. 
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dM/du = -T m = -\- 2 T uu (FNI) . 

We find that 

rh = -f uu /2X 2 (FNI) . (85) 

Alternatively, we can derive this relation using the Vaidya-like solution constructed in Ap- 
pendix A. To that end, we recall that the ingoing fluxes of both energy and charge, as well 



as the semiclassical correction term T uv , vanish at FNI (where R — > oo) - cf. Eqs. ( 55 , 
|79|) . The only semiclassical correction terms which survive at FNI are the Hawking energy 
outflux T uu and charge outflux q tU . We can therefore represent the solution near FNI by the 
(charged) outgoing Vaidya-like solution. This exact solution provides the desired relation 
between rh and the energy outflux at FNI. Identifying f uu with T^ ux) in Eq. flATOp below, 
we recover the relation (|35"D. 

Substituting the above expression for T UU (FNI) in Eq. (|35D, we find 



rh = -K[k 2 JX 2 + (gq/R + ) 2 ]/2 = -K[X 2 (m 2 - q 2 ) + g 2 q 2 ]/2R 2 + . (86) 

To verify the consistency of the above results for rh and q (and, more generally, for 
the fluxes of energy-momentum and charge), we can calculate the rate of change of R + 
in two different ways. First, since the evaporation is very slow (corresponding to a large, 
macroscopic, BH), we can view the geometry as quasi-static. At each "moment" u, R + (as 
viewed by a distant observer) can be estimated by the momentary values of m and q, via 
the standard, static-solution relation R + = m + (m 2 — q 2 ) 1 / 2 . Alternatively, we can apply 
the constraint equation R vv — R tV S tV + T vv = to the null generators of the EH, and in 
this way to analyze the rate of contraction of R + . Since the evolution is slow (and hence, 
on time scales short compared to the BH's evaporation time, the geometry only depends on 
x = v — u to the leading order), dR + /dv at the EH must coincide with R + = dR + (m, q)/du 
at FNI, i.e. with mi? + m + qR +tQ . In Appendix B we calculate these two quantities and 
show they are indeed the same: 
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R+, V (EE) = R + = -KA 2 [A 2 (m 2 - q 2 ) - g 2 q 2 ]/(2k + R 2 ) . (87) 

Finally, let us compare our result (|84D for the Hawking outflux at FNI to the standard 
result obtained by CGHS (in the uncharged case). Taking the limit q = in Eq. (j84|) , one 
finds 

f uu = KX 4 m 2 /R 2 + = KX 4 /A (FNI , q = 0) . (88) 

Transforming this result to the asymptotically-flat u, v coordinates defined in section II 
(which is also the gauge used by CGHS), we find 

f m = KX 2 /A = NX 2 /9Q (FNI , q = 0) . (89) 

This is just one half of the Hawking outflux in the CGHS model. This difference is because 
the quantum matter field used in the present model is fermionic, whereas that used in the 
CGHS model is bosonic [fLOR . 

B. Evolution of the charge-to-mass ratio 

Equations ([83]) and (|8~6D form a closed autonomous system, which allows us to analyze 
the evolution of the charge-to-mass ratio upon evaporation. One finds 

dm/dq = (R+/2q) [(k + /X 9 ) 2 + q 2 /R 2 + ] . (90) 

Since the right-hand side only depends on q and m through q/m, this equation admits 
solutions of the form 

m = cq , (91) 

where c is a positive constant [to be determined from an algebraic equation based on Eq. ([)(]) . 
as we shortly show]. We shall refer to a solution of this form as the linear solution. To analyze 
this solution, we rewrite Eq. (|90| ) as 

dm/dq = (R+/2q) [(n + /Xg) 2 - l] + m/q . (92) 
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This form makes it obvious that for any < g < A/2, there exists a linear solution of the 
form (|9T|), with c defined by the algebraic equation K + {m = cq) = Xg, i.e. 



-i 



g/X=l-(l + Vl - c- 2 J (93) 
[cf. Eq. (0)]. Explicitly we find 

c = (1 - [(1 - g/\y l - 1]T 1/2 (° < 9 < A/2) . (94) 
Next we analyze the stability of the linear solution (|9l| , |94] ). To that end, we define 



5 = m/q — c, and write the evolution equation for 5 in the form 

dS ^ - ™ = (R + /2q) [{n + /Xgf - 1] = r(5) . (95) 



<i In g <ig q 



Note that k + is an increasing function of 5 [cf. Eq. (|32])]. Since Ag = K+(m/q = c) = 
k_i_((5 = 0), the term in squared brackets is an increasing function of S which vanishes for 
5 = 0. Therefore, T has the same sign as 5 , which means that \5\ is an increasing function 
of q. This implies that upon evaporation (q decreases), \5\ decreases. Namely, the linear 
solution ([TI|, |94j) is stable. Moreover, since the only zero of T is at 5 = 0, the linear solution 
is in fact a global atractor for any < g < A/2, provided that initially q > 0. [For small 
5, we can linearize T by T = f3S, with some constant f3 = f3(g) > 0. We then find that 
5 oc q 13 a mP ] 



The dynamical system (|83| , |86|) has another, trivial, solution: 

q/m = 0. (96) 
For < g < A/2, this solution must be unstable: As we have just found, the linear solution 



([H],[)4]) is a global atractor in this range for any initial q > 0. However, for g > A/2, for 
which the above linear solution does not exist (c is not real), the trivial solution q = 
becomes stable. To verify this, we define in this case 5 = m/q > , and analyze Eq. (|95|). 
Since now n + /\g < 1 (for any q/m), T(S) is always negative, meaning that upon evaporation 
S increases monotonically. Furthermore, since R + /m > 1, the quantity din 5/ din q = Y/5 
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is bounded above by the strictly negative number 7/2, where 7 = A 2 /4g 2 — 1 < 0. This 
means that upon evaporation (lng — ► —00), 5 gets unboundedly large positive values. Once 
S becomes large, we can use the linear approximation Y/5 = 7 (obtained by approximating 
R + = 2m and k + = A 2 /2), which yields 5 oc q 1 . This implies m oc q x2 ^ g2 , namely, 

q/m oc m 4 ^ 2 " 1 {g > A/2) . (97) 

It should be pointed out that this linear analysis of solutions with q/m << 1 , and particu- 
larly the result (|97|) , applies to any g. It indicates the stability of the solution q/m = in 
the range g > A/2 , and its instability in the range g < A/2. Thus, for g > A/2, Eq. fl9~7D is 
realized as the late-time, stable, asymptotic behavior. For g < A/2, however, even if initially 
q/m << 1, upon evaporation it grows according to Eq. (^) until the linear approximation 
breaks (provided, of course, that initially q > 0). Subsequently q/m converges to a nonzero 
value c -1 , as was discussed above. 

We conclude that for g < A/2, the charge-to mass ratio converges to a nonzero value, 

L _> v/l-[(l-^/A)-i -l] 2 (0 < g < A/2) . (98) 

This value is independent of the initial values of q and m (though it only holds if q is initially 
nonzero). However, for g > A/2, the charge-to mass ratio decreases as a power of m, and 
eventually approaches zero, as described in Eq. fl97D . 

V. SUMMARY 

We presented here a dilatonic two-dimensional model of a charged black hole. On the 
classical level, our model yields static charged BHs, characterized by the two parameters m 
and q (representing the BH's mass and charge). These static BHs admit an inner horizon 
instead of a spacelike singularity. Then we studied the semiclassical effects (on the back- 
ground of the above classical, static, BH solution), using the method developed in Ref. . 
We derived explicit expressions for the fluxes of charge and energy-momentum as a function 
of the "radius" R, both outside and inside the BH. 
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By analyzing the outflux of energy-momentum and charge at future null infinity (and 
also the influx at the EH), we calculated the evaporation rate of both m and q, as a function 
of these two quantities. This yields a system of two coupled first-order differential equa- 
tions, i.e. q(m,q) and m(m,q) [Eqs. (|S3|) and fl86|), respectively]. We then analyzed the 
evolution of the ratio q/m upon evaporation. Depending on the value of the electromagnetic 
coupling constant g (recall g = VSgA), there are two different regimes: For g > A/2, upon 
evaporation q/m decays to zero as described in Eq. (|97| ) above. On the other hand, for 
g < A/2, the charge-to-mass ratio approaches a non-zero constant given in Eq. (f^). This 
constant depends on g but not on the initial values of m and q (provided that q is initially 
nonvanishing) . Note that this final charge-to-mass ratio approaches extremality (q/m = 1) 
at the limit g — > 0, and q/m — > at the limit g — > A/2. 

As was explained in the Introduction, no attempt was made here to investigated the final 
state of evaporation. The analysis throughout this paper was restricted to the macroscopic 
phase, i.e. to the stage where the mass is sufficiently large. This restriction is necessary for 
the validity of the adiabatic approximation: This approximation assumes that in evaluating 
the semiclassical fluxes (more specifically, when solving the field equations for the quantum 
variables Y and Z), m and q may be regarded as fixed parameters (and the background 
geometry may be approximated by the corresponding static BH solution). Clearly, this 
approximation is only valid as long as the change in m during a dynamical time scale is 
much smaller than m itself. The dynamical time scale (expressed in terms of u and/or v) 
is of order 1/k + , which is typically of order ~ A~ 2 . The mass evaporation rate rh is of 
order KX 2 (recall K = N/24). Thus, the macroscopic phase - the domain of validity of the 
adiabatic approximation - is given by 

m->K. (99) 
In this domain the dilaton field <fi outside the BH satisfies 

e~ 2(t> > K . 
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(This also holds inside the BH, in the region R > R_ - provided that q/m is not too small.) 
This is known to be the "weak-coupling" domain in large- N dilaton gravity. It should also 
be pointed out that the curvature singularity found in Ref. (in the uncharged case) occurs 
at e~ 2 ^ = 2K, which does not occur in the macroscopic domain considered here. 
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APPENDIX A: VAIDYA-LIKE SOLUTIONS 

In this Appendix we describe the construction of the Vaidya-like solution in double-null 
coordinates. For concreteness we shall consider here the outgoing solution, which is useful 
for describing the geometry of an evaporating BH (in the weak-field region), but the ingoing 
solution can be constructed in a completely analogous manner. 

Consider first the case of an uncharged null fluid (i.e. q is a fixed parameter). The 
solution is uniquely determined by the function m{u). In analogy with the construction of 
the static vacuum-like solutions in section II, we define 

H (R, u) = X 2 [R - 2m(u) + q 2 /R] . (Al) 

The function R(u, v) is now determined by the ordinary differential equation 

R, u = -H(R,u). (A2) 

This equation is to be integrated along the lines of constant v (with each ingoing null line 
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having its own "initial value" for R), and this integration produces the function R(u, v).f\ 
Then, S(u, v) is given by 

e s = R >v . (A3) 

Differentiating Eq. (|A2|) with respect to v, one recovers the field equation ( |2"3"D for R, i.e. 
R,uv = Fe s . Next, differentiating Eq. (|A3|) with respect to u, one obtains S >u = F, and a 
second differentiation with respect to v now yields the field equation for S. 

To discuss the null-fluids content of the Vaidya-like solutions (or, more generally, the 
radiative solutions), it is useful to re- write the constraint equations as 

— T uu = R w - R tU S, u + Tj« , (A4) 
= T VV = R iVV - R tV S, v + . (A5) 

The components Tuu UX ^ and T^ lux ^ then describe the energy-momentum carried by the 
outgoing and ingoing null fluids, respectively. (The vacuum-like solutions then correspond 
to Tuu UX ^ = Tvl lux ^ = 0.) Differentiating Eq. ( |A2j ) with respect to v yields R yV S )V = R >m . 
Namely, the outgoing Vaidya-like solution is characterized by a vanishing influx: 

= . (A6) 

To express the flux in the u direction in terms of m(u), we notice that the above result 
S, u = F together with Eq. (|A|) yield R tU S, u = -FH and R )UU = -FH - dH/du [where 
dH/du = (dH/dm)(dm/du) ], and hence 

T { J! UX) = R u S,u - Run = dH/du = -2\ 2 m, u . ( A7) 

3 This construction fixes the gauge for the outgoing coordinate u, but leaves the gauge for the 
ingoing coordinate v uncpecified. In this construction v enters as a parameter which parametrizes 
the one- parameter set of solutions to the differential equation ( |A2| ) for R. (For example, one can 
take v to be the value of R on some "initial" outgoing ray u = Uq.) 
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Note that Tuu UX ^ is independent of v - this is an important feature of the Vaidya-like solutions 
[valid for any F(R)]. 

Strictly speaking, the field equations in the form ( ]23|j2^) assume that q is a constant. 
One can, however, immediately generalize it by allowing q to be a function of u and v. When 
considering the (outgoing) Vaidya-like solutions, it is most natural to assume that q (like m) 
depends on u only. Physically, this would correspond to a model with an outflow of charged 
null fluid. This generalization is important because in our model the Hawking outflux is 
indeed charged. 

The generalization of the Vaidya-like solution to the charged null fluid case is straight- 
forward. One simply replaces Eq. (|A1|) by 

H(R, u) = X 2 [R - 2m(u) + q{uf / R] . (A8) 



The rest of the above construction, Eqs. (|A2|) and (|A3|) , are unchanged. The solution is 
uniquely determined by the two functions m(u) and q(u), which describe the outflux of mass 
and charge, respectively. One can verify that the evolution equations (0,|^) are satisfied, 
as well as Eq. ( |A6|) . However, the energy-momentum content of the outgoing flux is now 
modified: 

Tg lux ^ = R, U S, U - R, uu = dH/du = -2A 2 (m, u - q q jU /R) . (A9) 

Note that T$ ux) is no longer constant along lines of constant u . This has a simple physical 
interpretation: The Lorentz force acting on the charged outflux does a work on it, and 
changes its energy- momentum content. (This situation is well known in the context of the 



four-dimensional, spherically symmetric, charged Vaidya solution (TTJ; see the discussion in 



1 372 1 and ]T3|] .) Note, however, that at FNI the term q u /R vanishes, and one again obtains 

= -2\ 2 m tU (FNI) . (A10) 

We can verify this result by relating Tuu UX ^ to the rate of change of the Bondi mass M. In 
terms of the asymptotically-flat null coordinate u, these two quantities are related by 
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dM/du = -f fifi . 

Substituting M = 2Am, u = Xu, and T^u = X~ 2 T UU , we find 2m >n = dM/du = —\~ 2 T UU , 
which conforms with Eq. (|A10 ). 



The Vaidya-like solution may be interpreted as a slowly varying, quasi-static solution, 
which is described by the vacuum-like solution fl2~T|) - except that the BH's mass and charge 
are slowly evaporating. This interpretation is meaningful as long as the evaporating BH is 
macroscopic (i.e. the evaporation time scale is much larger than the dynamical time scale). 
Note that in this quasi-static limit the coordinate u used in the above construction of the 
Vaidya-like solution coincides with the Eddington-like coordinate u of the static, vacuum- 
like, solution ( |27D . That the gauges of these two solutions agree can be seen, for example, 
by recognizing that R yU = —H in both solutions. 

APPENDIX B: RATE OF CONTRACTION OF THE EVENT HORIZON 

In this Appendix we shall calculate the contraction rate of R + in two different ways, as 
outlined in section IV. The equality of the two results may serve as a consistency test for 
the expressions derived above for rh and q. 

First, since the evaporation is slow, we may assume that at each moment u, R + is given 

by 

R + (u) = m(u) + [m(uf - q{u) 2 ] 1/2 . 
Taking ^-derivatives of all quantities at FNI, one finds 

R+ = m [l + m/(m 2 - q 2 ) 1/2 ] - qq/{m 2 - q 2 ) l/2 . (Bl) 

Noting that 

1 + m /( m 2 _ g2) l/2 = R+/{m 2 _ g 2y/2 = X 2 /k+ 

[cf. Eq. (0)], we can re- write Eq. ( pl|) as 
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k + R + = A 2 (m — qq/R + ) . 

Substituting the above expressions for rh and q, one finds 

K + R + = -KX 2 [X 2 (m 2 - q 2 ) - g 2 q 2 ]/2R 2 + . (B2) 

In the second way, we apply the constraint equation (Q) to the null generators of the 
EH, and use it to analyze the rate of decrease of R + with v. Since the evolution is assumed 
to be slow, we can neglect the term R jVV and write 

S,v R+,v = T m (EH) . 



Using the background solution (|27|) one can easily verify that at the EH S jV = dH/ dR = 2k_ 
, and Eq. fl79|) yields 



T m = -K[k\ - (Xgq/R + ) 2 ] (EH) . (B3) 

Therefore, 

k + R +iV = f vv /2 = -K\ 2 [A 2 (m 2 -q 2 )- g 2 q 2 ] /2R 2 + (EH) . (B4) 

Comparing Eqs. ( |B2"|) and ( p4[) , we find that indeed R+, v at the EH is exactly the same as 
R + , as one should expect. 
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